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Abstract. We have studied systematically the free-electron laser in the context of high brightness pulsar radio 
emission. In this paper, we have numerically examined the case where a transverse electromagnetic wave is 
distorting the motion of a relativistic electron beam while travelling over one stellar radius (~ 10 km). For 
different sets of parameters, coherent emission is generated by bunches of beam electrons in the radio domain, with 
bandwidths of 3 GHz. Pulse power often reached 10 13 W, which corresponds with brightness temperature of 10 30 K. 
The duration of these pulses is of the order of nanoseconds. In the context of pulsar radio emission, our results 
indicate that the laser can produce elementary bursts of radiation which build up the observed microstructures 
of a few tens of microseconds duration. The process is sensitive mostly to the beam particles energy, number 
density and the background magnetic field, but much less so to the transverse wave parameters. We demonstrate 
that the operation of a free-electron laser with a transverse electromagnetic wiggler in the pulsar magnetosphere 
occurs preferably at moderate Lorentz factors 7 > 100, high beam density n > 0.1 7iGj(r*) where nGj(^*) is the 
Goldreich- Julian density at a stellar radius r», and finally, at large altitude where the background magnetic field 
is low B < W' 2 T. 
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1. Introduction 

Although radio pulsars have been discovered now for 
more than three decades, the generation of the radio 
signal which reveals their existence, is still a puzzle. It is 
believed that a pulsar is a rotating neutron star with a 
high magnetic field (B ~ 10 5 — 10 9 T). The radio pulses 
span a frequency range from ~ 100 MHz to ~ 30 GHz 
and have the remarkable property that, in a given pulsar, 
they vary dramatically from pulse to pulse, while the 
average pulse profile is extremely stable and unique for 
that pulsar. On one hand, the stability and uniqueness 
of the averaged pulse shape for every pulsar imply that 
the radio pulse is generated well inside the light cylinder 
(defined as the cylinder with radius R\ c = c/fi*, where f2» 
is the rotation frequency of the pulsar and c is the speed 
of light). On the other hand, the variability of successive 
pulses suggests that the radio emission process is strongly 
fluctuating and/or that the acceleration of particles 
responsible for the emission is highly non-steady. The 
observations indicate that the radiation is emitted more 
or less tangential to the open field lines from the magnetic 
poles and that it is polarised preferentially in the local 
plane of the magnetic field lines and perpendicular to 



it. Current belief is that radio emission is related to 
the development of a cascade of pairs of electrons and 
positrons on the open field lines anchored at the magnetic 
poles. 

An important property of the pulses is their 
high brightness temperature Tj,, given by Tb ~ 



10 34 K F Jy dl pc 



v, 



GHz ' 



where Fj y is the measured 



flux in Jansky at the frequency i^ghz for a pulsar at 
distance d^ pc and where r is the radius of the emission 
region (in meters). For characteristic values: F ~ 1 Jy, 
d ~ 1 kpc at v ~ 1 GHz, this implies that the brightness 
temperature ranges from 10 26 K to 10 30 K, depending on 
whether the radius of the emission region is taken to be 
the whole pulsar (r = ~ 10 4 m) or only the polar cap 



(r 



' pc 



tvCOt./c)- 1 / 2 ~ 10 2 m for P* = 2vr/0„ = Is). 
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If the emission process were incoherent this would imply 
the presence of energetic electrons (positrons) of indi- 
vidual energies £ ~ k B T/m c ~ 10 10 - 10 18 GeV, which 
are difficult to achieve in view of the maximum available 
voltage jump in a rotating magnetic star inside the light 
cylinder. Because of this, the mechanism responsible for 
the radio emission is assumed to be some form of coherent 
action. The emission can then be either an antenna or a 
maser process. If the radiation source is located on the 
open field lines above the polar cap a maser emission 
process is most likely powered by a high-energy elec- 
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tron/positron beam. Candidates f or the emission 



are maser curvature emissio n llLuo fc Melrose 
relativistic plasma emission l|Melrose fc Gedalir 



process 




norm al and anomalous Doppler inst ability lILvutikovl 
Il998l) . and linear acceleration emission l|Melroselll97S|) . 

Despite the efforts in the investigations over the past three 
decades, the mechanism responsible for the radio emission 
is still unknown. 

The process studied here is a specific type of the free- 
electron laser (FEL), a laboratory device which pro- 
duces coherent radiation. In a FEL, a beam of rela- 
tivistic electrons, with velocity v, passes through a pe- 
riodic, electromagnetic field (called wiggler) and radiates 
at the resonance frequency v rcs = w ros /27r = ±(2 r ) 2 (uj VJ — 
k w v))/2ir, where w w and fc w are the wave frequency and 
the wavenumber respectively associated with the wave- 
length A w of the wiggler. Under appropriate conditions the 
radiation will be enhanced by the particles that follow, due 
to the bunching of the par ticles. Previou sly mentioned lin- 
ear acceleration emission dMelroselll978i) and co herent in- 
verse Compton scattering ( Sc hopper et al.lE002l) are both 
forms of this mechanism. They differ from the case here in 
the sense that, in this study, particles will undergo a small 
transverse displacement, whereas in the other two cases, 
the distortion is along the path of the particles. Therefore, 
the solution proposed here is only applicable for the region 
where the background magnetic field is sufficient small, in 
comparison to the magnetic field amplitude of the wiggler, 

1. e. high up in the magnetosphere. 

Applying the FEL concept to the pulsar, we will inves- 
tigate as potential wigglers high-frequency Alfven waves 
and other potential periodic structures in the pulsar mag- 
netosphere. Such waves might be generated by a beam 
instability of the (remnant of the) primary beam in the 
ambient secondary pair plasma, or by the inhomogene- 
ity of the pair plasma in which faster particles run into 
dense clumps of pairs lIUsovl Il987t lAsseo k, Melikidzel 
Il998l IMelikidze etaP l2000|) . The efficiency of the FEL 
interaction between these waves and the beam of pri- 
mary/secondary particles is investigated by doing numer- 
ical simulations. 

In section 2 the basic concept of the free-electron laser 
mechanism is presented in detail before it is applied to 
the pulsar magnetosphere. A description of the simulation 
code is given is section 3, and a summary of free-electron 
laser parameters is found in section 4, In section 5, we 
will present the results of the numerical simulations of 
this process under pulsar magnetosphere conditions. The 
conclusions and a discussion are given in section 6. 
Throughout the paper, SI units and Cartesian coordi- 
nates, in which x is the position vector, with local mag- 
netic field parallel to z are used. 

2. The Model: A FEL in the pulsar magnetosphere 

In this section the basic theory of FEL operation is ex- 
plained. 



Note that the FEL mechanism that we treat here is an 
antenne process in which bunching occurs in space ("re- 
active" ) . 

2.1. Free-electron laser theory 

2.1.1. 1-particle trajectory and resonance frequency 

The relativistic equation of motion for a particle in an 
electromagnetic disturbance (wiggler), with electric field 
given by E w and magnetic field by B w , is given by: 



d'ymfBc 
dt 



= q[E w +/3cX B w ] 



(1) 



where the electromagnetic field of the wiggler is described 
by: 

E w (x, t) = E w cos(fc w z — iv w t)x 
B w (x, t) — B w cos(fc w z — oj w t)y 



(2) 



and (3 W denotes the wiggler phase velocity normalised to 
the velocity of light c. For a relativistic particle, with ini- 
tial Lorentz factor 70 = (1 — /3 z o) ^ 1, the first order 
approximation of Q, where we neglect the particle's en- 
ergy change, becomes: 

— - = uj w 1 - — — cos(k w z - w w t) (3) 

dt 70 V Pw J 



The solution of this equation is: 
V2K 



fi x = 77 sin(fc w z - uj w t) 



7 



(4) 



where j] = q/\q\ and the strength of the wiggler is re- 
expressed in a dimensionless parameter K, defined as: 



K = 



qB w X v 



(5) 



Thus, the effect of the wiggler on the particle orbit is 
adding a small sinusoidal motion in the direction of the 
wiggler 's electric field, where the periodicity of this mo- 
tion for an observer in the laboratory frame is given by 



A,-, 



A™ 



(PzO - /3w)c 



(6) 



Up to the second order, the (angular) frequency of the 
radiation emitted by a particle moving in this disturbance 
is derived from (u> — k • v) = ±(cj w — k w v z ): 



27 2 (w w - k w v z ) 



1 + K 2 



(7) 



where we have used u> - 
beamed forward and 01 



k • v 

= (3 2 - 



■0* = 1 



for radiation 

- 2 -K 2 j- 2 . 
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2.1.2. Collective effect: bunching 

The wiggler field causes a transverse motion of the par- 
ticles (0J). In presence of a radiation field E ra( j, B rac i oc 
sin(fcz — u>t), which is set up by particles radiating at the 
resonance frequency (JJJ , the particles that follow will move 
in a beat wave of the wiggler field and the radiation field. 
This beat is often called the pondcromotive wave, for its 
spatial energy density is highly non-uniform. The corre- 
sponding ponderomotive force is F p = q(3 x cx x B ra( j and 
is therefore proportional to sm((k + k w )z — (uj + uj w )t). The 
beat which is proportional to sin((fc — k w )z — {uj — u> w )t) is 
superluminal, and, therefore, has no immediate particle- 
wave resonances. 

This force acts on the beam in the z-direction and drives 
a longitudinal current 5J Z = qnv p , where q is the elemen- 
tary charge of the beam particles, n is the number density 
of the particles and v p is the velocity which is induced by 
the motion in the ponderomotive wave. As a result, like 
F p , 8J Z is proportional to sin((fc + k w )z — (to + tu w )t) as 
well. Related to 6 J z is a density modulation according to 
qdSn/dt — — V-<5J z z. This density modulation is observed 
as bunching of the particles. Note that the bunches occur 
at the ponderomotive wavelength 

A P = 2nk p 1 (8) 
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Fig. 1. Values for which w w , k w and 7 give rise to radio 
emission between 100 MHz (lower solid line) and 10 GHz 
(upper) according to the resonance condition JJJ with 
K = 10. On the vertical axis is plotted fc w w — cj w , be- 
cause we assume the particle's velocity to be larger than 
the wiggler's phase velocity. 



where k p = k w + k. Due to this bunching, the transverse 
current is in phase with the ambient field (5J± = q6nv w oc 
sin(fcz — u>t)), i.e. particles in the bunch radiate in phase 
and enhance the ambient radiation field. 



2.2. Background Magnetic Field 

The effect of a background magnetic field B = Bqz lim- 
its the particle's motion perpendicular to the field. The 
equations of motion for a particle now become: 



d7& 
dt 



dt 

di(3 z 
dt 



V2Kuj w (l — — - ) sin(/c v 

V Pw/ 



z - uj w t) + n Bo (3 y (9) 
(10) 

(11) 



2.3. Application in the pulsar magnetosphere 

Our aim is to apply the mechanism to the pulsar mag- 
netosphere, while making as few assumptions as possible. 
As is common to most pulsar models, we assume that a 
large electric field is set up along the open field lines above 
the polar caps due to the fast rotation of the magnetised 
neutron star. This electric field pulls out and accelerates 
electrons, which, in the presence of a strong magnetic field, 
are in the lowest Landau level and thus move along the 
magnetic field lines. At some altitude above the pulsar 
surface - between a few polar cap radii and a few stellar 
radii - a dense plasma of electron-positron pairs is pro- 
duced, either from curvature radiation of the primaries in 
the strong magnetic field or from their inverse Compton 
radiation. Very likely, the entire process is only stationary 
in an average sense but highly variable in space and time 
on small scales with inhomogeneous distributions of pair 
plasma intermingled with primary beams. 



where K is defined by equation © and £Ib = \q\B/m is 
the non-relativistic gyrofrequency for a particle (mass m, 
charge q) in a magnetic field with strength B. For zero 
background magnetic field, the equations reduce to QJ. 
In the presence of a uniform background magnetic field, 
the first order solution for the particle's velocity compo- 
nents (the x- and y-component) are coupled to each other 
through flso- The particle is no longer free to resonate 
with the wiggler, but is bounded by the magnetic field. 
This restriction limits the particle's ability to resonate 
and, therefore, also the FEL laser action. 



In our model, we investigate radiation from a mono- 
energetic beam of electrons in the presence of a subluminal 
transverse electromagnetic wiggler. 

In theory, radiation of any frequency can be obtained 
by tuning cj w , k w and 7. For radiation at frequencies 
u Tes = (x> ros /27r in the range of 10 8 and 10 10 Hz, and paral- 
lel wave vectors k w || v, the required Lorentz factors range 
from unity up to 10 7 (Fig.^). 

Near the stellar surface, particles move one-dimensionally 
along the magnetic field lines due to extremely fast syn- 
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chrotron losses in the strong magnetic field. The quantity 3.1. Code General Particle Tracer (GPT) 



Pj_/Bo is not invariant anymore when the timescale at 
which the magnetic field changes, is small compared to 
one gyration period in the background magnetic field, i.e.: 



B 2tt7 
dB w /dt <<C ShT 



and in our case, using J3J): 



B < 2tt< 



'1 



Ay 



(12) 



(13) 



So, if the last inequality is satisfied, a particle moving in 
the background magnetic field acquires momentum trans- 
verse to the magnetic field. 

Next to consider is the synchrotron loss time, which needs 
to be larger than the FEL timescale when the latter pro- 
cess causes the pulsar radio emission. The synchrotron loss 
time T1/2, the time within which the particle looses half of 
its initial energy E — 7771c 2 , is defined as T1/2 = E/(2P), 
where P = r ) 2 0jJl 2 B /(67re c) is the power of synchrotron 
radiation emitted by a charged particle moving in a back- 
ground magnetic field Bo and /3j_ is the component of the 
normalised velocity perpendicular to Bq. By using the up- 
per limit for Bq (|13[1 . we find for T\/2 



A 



1 3e mc 

Tl/2 ^ 2 2V (f3 z0 - /3 W )2 K*>y 



(14) 



where we used f3± w K/j from J2J|. 

On timescales much smaller than T1/2, a change of the 
particle energy due to synchrotron radiation can be ne- 
glected. 

In our calculations, we used parameters such that firstly, 
the resonance is in the radio regime |7|. secondly, the in- 
equality (|13fl is satisfied so that particles can (and do) ac- 
quire transverse momentum, and finally, the synchrotron 
emission can be neglected ()14|l. 

3. Numerical Simulations 

This section describes how the ingredients for the FEL in 
the pair plasma of a pulsar magnetosphere are represented 
in our numerical simulations. The Coulomb interactions 
between the particles in the beam are neglected, due to 
the large Lorentz factors. This is because, although the 
beam particles generate an electric field radially outward 
E r , due to their relativistic speeds, the generated magnetic 
field (in the azimuthal direction) reduces the electric field 
by a factor of 1 — /3 2 , where f3 is the particles' speed. This 
results in a reduced radial Lorentz force F r = qE r /^ 2 . 
Furthermore, the role of the pairplasma is mainly in pro- 
viding the wiggler. Therefore, we model the presence of 
the wiggler by equations J2J rather than generating it in 
a consistent way. 

First, we give a brief introduction to the code. Then, the 
parameters for each ingredient in the case of A^-particlcs 
simulations are given. 



We will investigate radiation from a relativistic electron 
beam travelling through a wiggler, and the formation of 
bunches is tested by doing numerical simulations with a 
code called General Particle Tracer (GPT) 1 . 
This code solves the equation of motion of each 
(macro)particle in the time domain numerically. For each 
macroparticle, labelled i, the differential equations 



dt 
dxj 
dt 



m,j c 



[E(x,t)+v<(f) x B(x,t)] 



2 R 2 



(15) 



7?/3 



are solved numerically with a fift h -order embedded Runge- 
Kutta ODE solver IPress et all l)l992j) . Here 7 = (1 - 
f3 2 )~ 1 / 2 is the Lorentz factor, (3 = v/c, and E(x, t), 
B(x, t) are the electromagnetic fields in which the particle 
is moving at position r and time t. Because a macropar- 
ticle represents N q particles of elementary mass 77i e and 
charge — e, the fraction q/m in the equation of motion is 
the same as in the case of single electrons. The electric field 
consist of two parts, E(x, t) — E w (x, t) + E rac i(x, i), where 
E w (x, t) is the wiggler field and E ra( j(x, t) is the radiation 
field, of which the development is studied in time; the same 
holds for B(x, t). The space charge effects of the particles 
in the beam are neglected. The initial conditions of the 
particles (i.e. number of charged particles, mass, charge, 
initial distributions in coordinate and velocity space) as 
well as the spatio-temporal behaviour (including initial 
amplitude and phase) of the wiggler, are user-specified. 

To find the radiation in the simulations, a set of differential 
equations for the radiated energy are solved, additional to 
the differential equations (|15|) . To derive this set of differ- 
ential equations, we first note that there are two essential 
differences between a FEL in the lab and our astrophys- 
ical application: firstly, there are no reflective mirrors in 
the neutron star magnetosphere. For a FEL to operate un- 
der such conditions, it should be a high-gain, single-pass 
process. Secondly, there are no side walls bounding the 
FEL cavity in the magnetosphere. This is accommodated 
in our pulsar study by using Gaussian modes of which the 
field decreases in a gaussian manner to zero away from the 
axis. 

The radiation electric field is splitted into a set of Gaussian 
modes (see Appendix A): 



E rad (x,i) = ^2Ej(x,t)x (16) 

j 

E 3 (x,t) = AjTjcoBfa) (17) 



1 see http:/ /www. pulsar.nl/ for a description of the code and 
other publications with this code 
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where 



Wj {z) 



k 3 r\z 



arctan 



z o,j 



(18) 



0,(19) 



where we have chosen cylindrical coordinates (z, = 
(x 2 + y 2 ) 1 / 2 ), and the polarisation is chosen to lie in the 
^-direction, in agreement with the polarisation of the wig- 
gler which is expected to generate radiation polarised in 
the x, z-plane. Note that, there should be a small compo- 
nent of electric field in the z-direction as well. But since 
this is negligible, as compared to the x-component of the 
Gaussian mode, its interaction with the particles is ne- 
glected in the calculations. 

Further, Wj(z) — wq^JI + z 2 / 'z^ j is the waist, w is the 

waist at z = 0, and z j — (l/2)kjWQ is the characteris- 
tic distance where the wave starts to diverge, Aj is the 
amplitude of this wave at wavevector kj, ojj = kjC is the 
corresponding angular wavefrequency, and ^ is an arbi- 
trary phase. 

Note that Gaussian modes propagate initially only in one 
direction and, therefore, they are suitable to describe nar- 
row beams of pulsar radio emission. Figure |2] shows the 
wavefront of a Gaussian mode (dotted line) ; the solid lines 
are level curves at e~ 1: ~ 2 '" of the maximum amplitude 
(left: near-field, right: far-field). As the figure shows, the 
wavefront of a Gaussian mode changes from planar to 
spherical. The asymptotic opening angle <f> of this wave 
is given by: 



tan ( 



A 



2ttzq 



(20) 



In the GPT, the user specifies the minimum and the max- 
imum frequency of the Gaussian modes and the num- 
ber of modes in this range. The frequency intervals are 
equally spaced (given by AN m ) and each is represented 
by a Gaussian mode with wavevector k. The spectrum in 
this range is represented by the sum of these modes (equa- 
tion l|17|0. When the radiation of the pulsar is assumed 
to be due to a FEL, then this set of Gaussian modes de- 
scribes the radiation properly. 

To continue the derivation of additional differential equa- 
tions, we use conservation of energy: the energy gain/loss 
per unit time of the ith macro-particle with charge Qi = 
N q qi due to the interaction with the jth electric field com- 
ponent is given by: nq^i -Ej, which is balanced by the rate 
of change of the electromagnetic fields. The total power 
radiated P = dW\i g h t /dt is obtained from: 



E., 



dWu, 



hi 



dt 



= 



(21) 



where the summation is over all N particles and N m 
Gaussian modes. 

Instead of calculating the change in amplitude and phase 
of each mode directly, GPT rewrites the amplitudes and 





Fig. 2. Wavefronts (dotted lines) and level curves at e~" 
of the maximum amplitude of the Gaussian mode (solid 
lines, n = 1,2,..); left: near-field, right: far-field, z is in 
units of zo and rj_ in units of wq . 



phases into two other independent variables rrij and ny. 

. \ , " 

y j J' 

(22) 



V J = ^/m 2 +n 2 /AN n 
4>j = arctan(nj/mj) 

Inserting equations (1221) into the conservation of energy 
gives two differential equations for mj and rij (j = 1.. 
N m ): 

drrij 
~df ' 



\- N q q t AN m 



nWyeoL 



drij x - N q qiAN ni 

= 2^ — —Vx.iTi cos( 



TTWqEoL 



(23) 



where L is the cavity length. 

A factor AN m enters Eqs. (1221) and H'23|) because one mode 
represents a frequency interval. This gives the correct cal- 
culation oirrij and rij, which is dependent on AN m , and 



6 



P.K. Fung and J. Kuijpers: A FEL in the Pulsar Magnetosphere 



ensures that the resulting spectrum, represented by Aj's, 

is independent on the number of modes used. 

The wiggler part of the electromagnetic field is assumed 

to have a time-independent amplitude during passage of 

the beam, and is not included here, for we assume that the 

change in the wigglers' energy is small on the considered 

timescales. 

Summarising, the differential equations solved at time t 
by GPT consist of the 2iV equations: 

*A = -^[E(x,i)+v 2 (t)xB(x,i)] 
at rrijC 



dxj 
~dt 



-fit* 



1 



for the N charged particles' phase space positions, and the 
set of 2N m equations: 



N qqi AN r , 

^ 1TWq£qL 

N q q t AN r , 



■V x ;T k C0S( 



drrij 

~dT : 

—7T~ = V 2 T v Xti T k COS(9j) 

dt t-r 1 ttwqCqL 



for the total number of modes N m included in the sim- 
ulation, to calculate the amplitude and phase at each 
wavenumber in the specified range during the interaction 
with the particles. As a result we obtain the spectrum of 
the radiation in time. 

To retrieve the pulse in the time domain, both the elec- 
tric field and the power are Fourier transformed from the 
frequency to the time domain, and Parseval's theorem is 
applied to the results to ensure a correct transformation 
between time and frequency domain. 
Note that the results obtained with this code are inde- 
pendent of the number of time outputs, the number of 
macroparticles and the number of modes. 



4. Physical Parameters for Numerical Simulations 

Except in the second last run of the simulations, we 
used a beam of 200 macroparticles with a total charge of 
— 6 • 10 ~ 4 C. These are uniformly distributed in a cylinder 
with radius of one metre and length of three ponderomo- 
tive wavelengths JHJ, which varies from run to run. The 
corresponding electron number density is 1.59 • 10 16 m~ 3 . 
This is 0.23 times the Goldreich- Julian density at the 
stellar surface, given by nQjfV*) = 2eol~i • B»/q, for a 
reference pulsar with Q, = 2tt/P, where P = 1 s and 
B„ = 10 8 T. Except in the last run, each macro-particle's 
initial Lorentz factor is 7 = 1000. 

The wiggler is an electromagnetic disturbance propagat- 
ing in the z-direction. In the pulsar magnetosphere, the 
scale at which the physical parameters change, e.g. back- 
ground magnetic field, is of the order of z/r* (> 10 km). 



Therefore, the interaction between the beam and the wig- 
gler is simulated until the beam particles reach this dis- 
tance, so i cn d ~ 3.3T0 _5 s. The parameters for the wiggler 
are given in Table ^ (all labelled with "w" ) . The values 
for the magnetic induction and angular frequency are ob- 
tained from (J2J. 

In the first run, the wiggler's phase speed is 0.9c. The am- 
plitude of the magnetic induction B w is 3 • 10~ 3 T. The 
wavelength A w is 50 metres. These values are chosen such 
that the resonance is at radio frequencies, and the strength 
of the magnetic induction is chosen such that the induced 
transverse velocity is small compared to the initial axial 
velocity (if/7 < 1 in Equation @). 

In the subsequent runs we investigate the effects of varying 
each of these parameters as compared to the results of the 
first run. In run 2, the wiggler strength K is increased from 
five to twenty, in steps of K — 5 (correspondingly the mag- 
netic induction is increased in steps of B w = 1.5 ■ 10 -3 T). 
In the third run, the wavelength of the wiggler takes the 
values 25, 40, 50 and 100 m. In run 4, the velocity of the 
wave is changed from 0.2 to 0.9. 

In run 5, we study the effect of a non-zero background 
magnetic field. We took the most simple case is which the 
field in uniform (in z) and has a background magnetic field 
strength of B = 10~ 3 , 10" 2 , 2.5 • 10~ 2 , 5 • 10~ 2 T. 

In run 6 as we change the number density n of the elec- 
tron beam. In terms of the Goldreich- Julian density, the 
fractional density ranges from 0.4 to 0.1. 

Finally, in run 7, we lower the Lorentz factor of the beam 
particles (the resonance frequency depends on this as 
uj ics cx 7 2 ). The Lorentz factors are chosen: 7 = 1000, 
500, 250, 100. 



5. Results 

Figs. ISaTRun 1) toEfilRun 7) and Table □ show the end 
results of all runs, i.e. at t = 3.3 • 10~ 5 s. Rather than dis- 
cussing the runs sequentially, we present the systematic 
trends of our computations, using the figures as illustra- 
tions. Run 5, where a uniform background field is included, 
is presented separately. 



5.1. Bunching 

We found bunching of particles during their interaction 
with the wiggler in a number of cases. Run 1 clearly shows 
this bunching (Fig. I3a[) . Here, the particles are plotted 
in (z,ir)-projection at different stages of their interac- 
tion: bunch formation (t = 1.33 ■ 10~ 5 s), 'de-bunching' 
(t = 2.0-10" 5 s) and 're-bunching' (t = 2.67-10" 5 s). Note 
that the distance between bunches is A p as expected. 
An alternative illustration of the formation of bunches is 
found in Fig. I3dl where the average Lorentz factor per 
macro-particle 7 is plotted. Evidently, the formation of 
bunches corresponds to a steep drop of 7 (e.g. Fig. |H3| 
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Table 1. Input parameters and results for all simulation runs. Input parameters are: electron initial Lorentz factor 70, 
number density in the beam n in terms of Goldreich- Julian density at the pulsar surface ncj*, dimensionless wiggler 
strength K, wiggler wavelength A w , dimensionless wiggler phase speed /3 W and (uniform) background magnetic field 
Bo- Results from the numerical simulations are: central frequency i^cn, FWHM bandwidth of the radiation Av, peak 
power of the light pulse P m ax, and brightness temperature corresponding to the peak power Tb- The last four entries 
are the energy density of each ingredient in the simulation: beam eb, wiggler e w , background magnetic field ebo and 
pulse e p . Note that every time "-" appears, it means the entry has the value as in Run 1. In Run 5, n appears for 
u reB , Av and Tb, which means, that for those cases, the spectra are flat, no peaks are seen, and also the brightness 
temperatures which depend on the bandwidths are not calculated. 



the first bunch formation starts from t ~ 1.0 -10 5 s to 
1.8 ■ 10~ 5 s, and the second formation from 2.4 • 10~ 5 s to 
2.8- l(T 5 s). 

We notice that both the bunch duration and the beam 
energy loss is larger for the first bunch formation. 
Saturation is reached due to an increased velocity spread. 
Therefore, for the second bunch formation, the beam 
starts with a certain spread in velocity, and it reaches 
saturation quicker than in the initial bunching. 
For the beam number density of n = 0.23 «gj*j wc nn d 
a beam energy loss to radiation of ~ 5%, after the first 
time bunching. This number is only affected by n. The 
number of particles in each bunch scales with n, as does 
A7 (Run 6, Fig. I8cjl . which agrees with coherent losses. 
As to the starting time of the first bunching, we find 
the following: for the same beam and different wiggler 
parameters (/3 W and A w ), about ten oscillations of the 
particles in the wave as seen by the observer are needed 
before they start to bunch (Run 3, Figs, loci and Run 4, 
Fig. I6c|) . For the same wiggler, increasing the beam 
density or decreasing the Lorentz factor give rise to 
earlier occurence of bunching (Run 6, Figs. I8c| and Run 7, 
Fig. |ScJ. This suggests that the bunching sets in when 
the (incoherent) radiation reaches a certain level, because 
the ponderomotive force is dependent on the radiation 
field. 



As for the wiggler strength K, it determines the maximum 
beam energy loss during bunching. In Fig. !4cf Run 2) are 
shown (7, t)-plots for K = 5, 10, 15, 20. We have also run 
cases for K < 1. The maximum energy loss occurs when 
K = 1. 



5.2. Coherent radiation 



That the radiation is coherent is demonstrated by the 
power of the radiation pulse (Fig. I3b(l . At the top is pre- 
sented the simulation run for N = 200 macro-particles 
and at the bottom for TV = 2 macro-particles. The max- 
imum pulse power for N = 200 is P 2 oo = 9.09 • 10 13 W 
and for N = 2, P 2 = 9.75 • 10 9 W. If the radiation were 
incoherent emission from 200 macro-particles, we would 
expect P oc N, where N is the number of particles and 
(-P20o/-P2)inc = 100. However, the simulation shows that 
(-P200/-P2) — 9 ■ 10 3 , and thus P oc iV 2 , demonstrating 
that the radiation is coherent. 

Note that only the maximum of the first peak satisfies 
P200 oc iV 2 . This is not the case for the second peak in 
Fig. I3bl where now P200 ~ 2 • 10 13 W, which is less than 
the coherent case but much more than if it were incoher- 
ent. 
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5.2.1. Spectrum 

The radiation frequency in our simulations agrees well 
with the expected resonance condition (compare eq. 
and Table Therefore, by construction, the frequency of 
the emission is in the radio regime, between 1 ~ 10 GHz. 
In some cases, the central frequency is lower than the res- 
onance condition (Run 2-4: Fig. I4al f5al I6a|) . The shift in 
these cases is caused by multiple bunching where the av- 
erage Lorentz factor decreases after each bunching, and 
therefore, also the frequency as determined by the res- 
onance condition. As a result of multiple bunchings, the 
total bandwidth increases (Runs 3, 4, 6 in Figs. l5all6all8a|> . 
whereas the relative bandwidth Av /i> decreases. 

5.2.2. Pulse shape, power, brightness temperature and 
radiation cone 

The pulse shapes of the coherent radiation can be found 
in Figs. l3hl l4bl l5hl Ifihl l8bl and l9hl 

The number of temporal peaks corresponds to the num- 
ber of bunching times (not the number of bunches that 
is formed); e.g. in Run 1 (N = 200) the beam particles 
bunched twice (Fig. I3d|) . and the pulse has two peaks 
(Fig. I3b(l . The maximum pulse power always occurs at 
the first time bunching. The following bunch formations 
result in less powerful pulses (though not completely in- 
coherent). The maxima of the pulses have the same order 
of magnitude for most runs (Tabled, except for 7 = 100, 
where the coherent emission becomes less efficient (and 
also for Bo = 0.05 T, but this is discussed in the next 
section). The optimum efficiency of the emission power 
is reached for 7 ~ 300 as follows from the ratio e p /eb 
(Table □ Run 7: Fig. Hp 

The structure within each pulse reflects the wiggling mo- 
tion of the particles (e.g. Fig. |HE1 where A w = 100 m). 
The typical duration of the pulse is a few nanoseconds. 
The brightness temperature Tb is derived using 



I = 



W/At 
AAvAQ 



(24) 



where / is the radiation intensity, fee is Boltzmann's con- 
stant, W/At is the power, which we take as P ma x, A is 
the emitting surface, which corresponds to the beam cross- 
section, Av is the bandwidth, and AO. = ir tan 2 <p the solid 
angle into which the radiaton is emitted H20(l . 
Converting the maximum power into the brightness tem- 
perature, we find that, again, for most cases the brightness 
temperature is approximately 10 30 K (Table QJ. 



5.3. Background magnetic field 

An ambient magnetic field in the z-direction (Run 5) 
can have a drastic effect on the bunching (Fig. [7c|l . For 
a background field just smaller than the wiggler field, 
B = IQ- 3 T, the particles' behaviour shows no difference 
with the simulations where the background field is absent. 
For Bq = 10~ 2 T, we still observe bunching. For larger 



values of Bq, however, bunching disappears completely. 
Actually, such behaviour is expected from equation (|13|) : 
since the motion of the particles is no longer free in the 
transverse direction , the ponderomotive force is less ef- 
fective and the resulting bunching is less pronounced. For 
all runs, the angle between the particles' velocity and the 
total magnetic field B = Bq + B w , which is proportional 
to v ■ B is indeed found to be nonzero. 
Since only the cases where Bo < 10 -2 T show bunching, 
these are discussed further. The spectrum for this run is 
plotted in Fig. [7a] Only for B = 10~ 3 T, the charac- 
teristics of the spectrum are similar to the cases where 
Bo = T. For B = 10~ 2 T, the resonance frequency 
shifts to v res = 2.18 GHz (as compared to the expected 
value of 11.9 GHz). The bandwidth is 0.61 GHz, and is 
much smaller than before. 

6. Discussion 

We have investigated a specific form of a single-pass free- 
electron laser process as a possible mechanism to produce 
high brightness radio emission of pulsars. We have investi- 
gated the operation of a FEL in the presence of a wiggler 
which consists of a transverse electromagnetic disturbance 
as it is being overtaken by a relativistic electron beam. 

We have shown that, in principle, a FEL can operate in 
the pulsar outer magnetosphere in the presence of a trans- 
verse wiggler. The deciding factors for particles to bunch 
and emission to be coherent within 3.3 ■ 10~ 5 s, are the 
following: a large beam particle density n > 0.1 ncj^*), a 
Lorentz factor of the beam particles 7 > 100 for K — 10, 
and a small background magnetic field Bo < 10~ 2 T. The 
brightness temperature of the pulse depends sensitively on 
these parameters (Table H). 

The required beam particle density together with the wig- 
gler's parameters mainly determine the timescale on which 
particles start to bunch. Bunching occurs after about 10 
times the transit time through the wiggler (i.e. 10 t e g, 
eq. Ell. This seem to set the level of inchorent radiation, 
which is then large enough for the ponderomotive force to 
act on the beam particles. Then the timescale over which 
a bunch stays together is about one i e ff, i-e. coherent ra- 
diation only occurs during this period. Due to an increase 
in axial velocity spread, debunching occurs. Although par- 
ticles rebunch, the associated pulse is weaker due to the 
velocity spread in the bunch. In the pulse profile, the first 
bunching results in one pulse. Then for larger times, parti- 
cles bunch more often and the pulse acquires more peaks, 
which are less powerfull. 

The overall fractional energy losses of the beam are ~ 5%. 
Together with the travel time, the calculations demon- 
strate that when applied to the pulsar magnetosphere for 
B < 10~ 2 T, the FEL interaction can produce coherent 
radiation. For a dipole model, this implies R 3> 2 • 10 3 i?*, 
where we used -B* = 10 8 T, i.e. in the outer magneto- 
sphere. 

While the particles are bunched, most radiation is at ^ rcs , 
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which in our runs varies from 0.7 GHz to 90 GHz as in 
the pulsar radio emission. The lower frequency is obtained 
with 7 = 250. 

The particles lose their energy most efficiently for K = 1. 
The measured FWHM bandwidth for the radiation Av 
ranges from 0.53 GHz to 14.5 GHz, and is broadband 
Avjv ~ 0.07 — 0.5. The width of the spectrum highly 
depends on the number of bunching events occurring dur- 
ing the simulation. After the first bunching, the band- 
width is of the order of 1.5 GHz; After the second -less 
effective- bunching, Av as 2.9 GHz and so on. Of course, 
the central radiation frequency agrees with radio pulsar 
emission by construction. However, it is interesting that 
the bandwidth is found to agree with pulsar radio obser- 
vations where average pulses and microstructures are ob- 
served from about 100 MHz to more th a n 10 GHz (average 
pulses : e.g.lLvne fe Manchester! lll988h: iManchester et al.1 
(|l996(k lLvneetal.1 lll99Sl): iD'Amico et all (Il998h. mi- 
crostr uctures: e.g. Ifiickett et al.1 l(l975h : iBartel fc Sieberl 

As for the brightness temperature Tb, which is a big ob- 
stacle for most radiation processes, we found a brightness 
temperature Tb at pulse maximum of 10 30 K for /3 W = 0.9 
and Tb = 10 29 K for a wiggler phase velocity (3 W < 0.8, 
similar to the observationally derived pulsar brightness 
temperatures. 

Apparently, the FEL is able to produce the required high 
brightness pulse. Further, the characteristic opening an- 
gle, given by tan0 w \/(2itwq), is <f) < 1°. Again, this 
agrees with the observed values in pulsar radio emission, 
as estimated from the microstructure duration relative to 
the pulsar rotation period: ~ 100^s/P x 360° < 1° for 
P = 0.1s. 

The radiation pulse has a duration of 2 ns. The short- 
est elements of radio emission measured from pulsars 
are microstructures. These are quasi-periodic structures 



polarized llRickett et al.1 Il975 


ICordes & Hankinsl 


lLanee et alJll998HPoDov et al. 


|2002l. 



Since the starting time for the particles to bunch is 10i c ff = 
10A W / (Pz — Av)j the faster the wiggler wave is, the longer 
it takes for the particles to get bunched. Also, the charac- 
teristic pulse duration t e g increases. This is clearly shown 
in Run 4, where we varied the phase velocity of the elec- 
tromagnetic disturbance between 0.2c and 0.9c. The com- 
putation time limitation forced us to consider only such 
relatively low values for electromagnetic wigglers in the 
relativistic outflow from a pulsar magnetosphere, but one 
can see from the particles' behaviour in Figure I6cl that 
even more relativistic electromagnetic disturbances would 
lead to longer timescales, e.g. for j w = (1 — /5 2 / ) _1/ ' 2 = 
100, the timescale would go up by a factor 5000, and be- 
come comparable to the observed micropulse durations. 
(The same argument holds for a lower beam number den- 
sity n, Fig.|Sc| i.e. when n < O.lncj*, bunching occur at 
t > 3.3- 10~ 5 s). 

Finally, we note that the coherent emission from electrons 



and positrons in a transverse wiggler add constructively. 
Therefore, the total number of electrons and positrons in 
a bunch, and not the charge excess, determines the emis- 
sion. 

In summary, the operation of a single-pass high-gain FEL 
with a transverse electromagnetic wiggler within the pul- 
sar magnetosphere in the radio regime, requires a mono- 
energetic beam of electrons/positrons at moderate Lorentz 
factors in a sufficiently small background magnetic field. 

We expect a FEL process to be also possible much 
nearer to the pulsar, where the magnetic field strength 
is large and the dynamics of the particles are one- 
dimensional, when longitudinal instead of transverse wig- 
glers are used. This would be the domain of single-pass, 
high-gain Cerenkov FELs, and a next logical step to study. 
Strong Langmuir turbulence h as been studied in this con- 
text by (|Schonper et al.ll2002|) . By using a Particle- in-Cell 
method, these authors show that electron scattering on 
Langmuir turbulence, which is excited in a self-consistent 
way, results in high power output of radiation. 

Acknowledgement: This work has been funded un- 
der a joint research project between the Centre for 
Plasma Physics and Radiation Technology (CPS) and the 
Netherlands Research School for Astronomy (NOVA) . We 
want to thank the anonymous referee for his/her com- 
ments. JK likes to thank Dr. Marnix van der Wiel and 
Dr. Theo Schep for stimulating interest in the project. P.- 
K. F likes to thank Bas van de Geer, Marieke de Loos for 
their support in using the code GPT and Cyrille Thomas 
for the description of the Gaussian modes. Also, we want 
to thank Kees van der Geer for setting up the simulations. 



Appendix A: Derivation of Gaussian modes 

T he full derivation of the Gaus sian modes can b e foun d 
in lMilonni fc Eberlvl l)l988ft and lErikson fc Singhl {1994 . 
As the radiation is produced by a narrow beam of rela- 
tivistic particles the radiation propagates mainly in one 
direction (z-axis) . As to the dependence of the radiation 
in the (x,;y)- plane, we assume that, similar to a laser beam, 
the intensity has a Gaussian form inside the FEL; i.e. 
I± cx e~( r - L l w °> where r±_ = yj x 2 + y 2 and wq denotes 
the transverse distance where the intensity drops to 1/e 
of the peak value at z — (wq is called the waist). 
The electromagnetic field which has these properties and 
satisfies Maxwell's equations is the Gaussian mode. The 
derivation of Gaussian modes starts from Maxwell's equa- 
tions in vacuum, which result in the (vector) wave equa- 
tion: 



V 2 E(r,i) 



1 <9 2 E(r,<) 
<? dt 2 



= 



(A.I) 



For vacuum waves (uj 2 = k 2 c 2 ) propagating mainly along 
the z-direction, the paraxial approximations can be ap- 
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plied: 

<9 2 E(r) 

dz 2 
d 2 E(r) 



< 2k 



<9E(r) 



dz 2 



dz 

< 2fc|E(r)| 



(A.2) 
(A.3) 



With these, the vector wave equation becomes the vector 
paraxial wave equation: 

<9E(r) 



VxE(r) + 2ik- 



dz 







(A.4) 



where 



d 2 /dx 2 + d 2 /dy 2 . This equation is valid for 
each component of E(r) = E x (r)x + E y (r)y + E z (r)z. The 
same holds for B(r). 

The general solution for each component of this equation, 
which is axisymmctric, is the following: 

ikr\ 



u{r± , z) = a exp(— P(z)) exp(- 



;) 



(A.5) 



2Q{z)' 

where a is an amplitude and P(z) and Q{z) are complex 
functions which specify the longitudinal and transverse 
mode behaviour. These can be retrieved by putting (|A.5|I 
into l)A.4p . which results in two differential equations for 
P(z) and Q(z): 

dQ 

dz 

£ = i (A.7) 



= 1 



(A.6) 



The solutions are simply given by: 

Q(z) = z + q 

P(z) = -i\n(z + q ) 



(A.8) 
(A.9) 



where qo is an integration constant. Recall that function 
Q(z) gives the transverse behaviour of the mode and that 
the solution is paraxial. Therefore, Q{z) can be written in 
terms of a radius of curvature R(z) and a width w(z): 

1 1 I -2i 

~Q{z) ~ ~z~+~qo ~ ~R{z) + Tm^Jz) ( A °> 

Furthermore, we assume that at a reference point z = 0, 
the mode wavefront curvature is i?(0) = oo, so that: 

ikw 2 (0) 



<7o = 



Together with: 

P(z) = -iln(z + iz ) = -iln[(z 2 + z 2 )e i0 ] 
= —i ln[z 2 + z 2 ] + arctan(zo/z) 



zq = 
tan <h = 



kwQ 



the general solution becomes: 
u(r±, z) = A — r^x 



w{z) 



Re 



-ikz-\-iu)t—ik - 



-ief> 



(A.ll) 

(A.12) 
(A.13) 

(A.14) 
(A.15) 



(A.16) 



where 

w{z) = 
R(z) = 




(A.17) 

(A.18) 
(A.19) 



where r± and w are as defined above, w(z) is the waist 
at z and A, k, uj are the amplitude, the z-component of 
the wave vector and the frequency of the electromagnetic 
wave respectively. The curvature radius of the wavefront is 
given by R(z) = z+z 2 /z, where z — (ttwq)/X is (roughly) 
the separation between the near and the far field. The 
wavefronts of this mode change from planar in the near 
field to spherical in the far field (see Fig.EJ). 

The opening angle ip of the intensity of the mode is given 
by: 

-w(z) _ 



tan( 



1 



V2- 



X 
2ttzq 



(A.20) 



Within distance of length zq, wavefronts of the Gaussian 
modes are considered planar. Different wavelengths and 
different zq result in different opening angles. In table lA.il 
is listed the range of opening angles when we consider 
radio waves with wavelengths between 3 cm to 30 m and 
some typical distances in the pulsar magnetosphere. 



z <f> 


polar cap radius R pc = 


100 m 


24' - 3°57' 


stellar radius R„ = 10 4 


m 


2'23" - 24' 


light cilinder i?i c = 10 7 


m 


4". 5 - 45" 



Table A.l. The range of opening angles <fi for different 
Zq and wavelengths between 3 cm and 30 m according to 



equation 
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Fig. 3a. Run 1: (z,x)-behaviour of the particles during 
their passage in the cavity (from left to rigth; from top to 
bottom). On the horizontal axis is plotted the z-coordinate 
of all particles, subsequently for z = 2, 4, 6, 8 and 10 km. 
Except for t = 3.34- 10~ 5 s, the horizontal spacing between 
the long ticks is 0.02 m (e.g. Az w 0.08 m for t = s). In 
the last plot, the spacing between the long ticks is 0.05 m. 
Bunching clearly occurs at t = 1.33- 10~ 5 s, 'de-bunching' 
at t = 2.0 • 10~ 5 s and 're-bunching' at t = 2.67 • 10~ 5 s. 
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Fig. 3b. Run 1: Power of the pulse at the end of the 
simulation t = 3.3 • 10~ 5 s (top: number of macroparti- 
cles N=200; bottom: N=2). On the horizontal axis is also 
plotted t, this is not the simulation time, but the dura- 
tion of the pulse as a distant observer would measure. 
Negative t corresponds to earlier arrival. The maximum 
power for = 2 (bottom) is P 2 ,ma X = 9.75 • 10 9 W. Since 
P 2 oo « N 2 P 2 , this implies the radiation pulse at the top 
is coherent. 
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Fig. 3c. Run 1: Spectrum at t = 3.3 • 10~ 5 s. The cen- 
tral frequency is at v = 11 GHz, with FWHM of Av = 
2.42 GHz. The peak on the right, centered on v Tea = 
11.9GHz, is caused by the first bunching at t = 1.33-10 s. 
And the peak on the left corresponds to the second buch- 
ing, where 7 dropped 5%, resulting in a lower v Ies . 




Fig. 3d. Run 1: Time development of average Lorentz fac- 
tor (i.e. per macro-particle) in two simulations, one with 
N = 2 (top) and the other with N = 200 (bottom). The 
N = 2 case reflects the behaviour of particles radiating 
incoherently, with constant energy loss over time as com- 
pared with the bunching case, where there is a steep drop 
in energy (A7 = 50) of the beam particles between 1.0 
and 1.8 TO" 3 s. Between 2.4 • 10~ 5 s and 2.8 • 10" 5 s an- 
other drop in beam energy occurs, which corresponds to 
a second bunching of the particles. This time, the beam 
energy loss is less than for the first time, and the duration 
is much shorter. The pulse power (Fig. I3bj) reflects these 
characterics where the first pulse is more powerful and has 
a longer duration than the second pulse. 
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Fig. 4a. Spectra for Run 2 plotted over each other; we 
have K = 5 (solid), 10 (dotted), 15 (short-dashed) and 20 
(long-dashed). The resonance frequency shifts according 
to Q.The bandwidth increases with increasing v res . 
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Fig. 4b. Pulse power for Run 2 plotted for different K. 
The maximum scales roughly as K~ 2 . 
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Fig. 4c. Run 2: (7,t)-plot for K — 5, 10, 15, 20 (labelling Fig. 5a. Run 3: Spectrum at the end of the simulation for 
as in Fig. I4a|) . Obviously, for all K, the beam particles different A w = 100 m (solid), 50 m (dotted), 40 m (short- 
bunch twice (also, see Fig.l4b|l. corresponding to the max- dashed), 25 m (long-dashed). Spectral broadening is due 
imum power, the beam energy loss is largest here for to more bunchings occuring during the simulation. 
I\ 5. 
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Fig. 5b. Pulse power for Run 3 plotted for different A w 
values. Maximum powers are of about the same order of 
magnitude. Only the number of pulses differs for each run, 
due to number of bunching times. 
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Fig. 5c. Run 3: The average Lorentz factor of the Fig. 6a. Run 4: Spectrum for (3 W — 0.2 (solid), 0.5 (dot- 
macroparticles during their passage in the cavity for A w = ted), 0.8 (short-dashed), 0.9 (long-dashed). 
100, 50, 40, 25 m (labelling as in Fig.lSaj). The starting 
time of the first bunching for each run is about 10 t c ff ©• 
The energy beam loss after the first bunching is about 5%. 



i ! i i i i ! i i i i 

L 


1 Q 


1 1 1 1 1 1 
-i . . . i - 


M Ml II Ml II 

£ w =0.2 ~ 


1 1 1 M 1 1 

\ 


1 | 1 1 1 1 | 1 1 1 1 

W =O.5 ~ 



-1 1 -22 
t (10- 9 s) 



-1 1 
t (10- 9 s) 



Fig. 6b. Pulse power for Run 4 plotted for different /?„ 
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Fig. 6c. Run 4: The average Lorentz factor versus the av- 
erage distance z of the particles in simulations with chang- 
ing /3 W = 0.2 , 0.5, 0.8, 0.9 (labelling same as Fig.lfiajl. The 
starting time of the first bunching is 10i c ff- 



Fig. 7a. Run 5: Spectra for (from top to bottom along 
the left vertical axis) B = 10~ 3 T (solid), 10~ 2 T (dot- 
ted), 0.0025 T (short-dashed), 0.050 T (long-dashed). For 
Bo = 10~ 3 T, the spectrum shows the same properties 
(central frequency and bandwidth) as in the case where 
the background magnetic field is absent. For Bq = 10 -2 T, 
the central frequency shifts to 2 GHz, whereas the band- 
width drops to 0.61 GHz. 
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Fig. 7b. Pulse power (logarithmic) for Run 5. The power 
of the radiation becomes smaller for Bq = 0.025 and 0.05 
T. 
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Fig. 7c. Run 5: The average Lorentz factor versus the av- 
erage distance z of the particles in simulations with back- 
ground magnetic field for B = 10~ 3 , 1CT 2 , 2.5 • 1CT 2 , 5.0 • 
10 -2 T (labelling as in Fig. I7a(l . The last case coincides 
practically with the horizontal line 7 = 1000. Only for 
B = 10~ 3 T and B = 10~ 2 T the beam particles show 
FEL action. 



Fig. 8a. Run 6: Spectrum (from top to bottom) for 
n/n,G3* = 0.1 (solid), 0.2 (dotted), 0.3 (short-dashed) and 
0.4 (long-dashed). The spectrum broadens as the number 
density increases, as expected from Fig. 02J which shows 
single bunching for n/riQj — 0.1, but multiple bunching 
for n/riQj > 0.1. 
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Fig. 8b. Pulse power for Run 6 plotted for different num- 
ber densities relative to ncj*. As expected, the maximum 
power scales as n 2 . 
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Fig. 8c. Run 6: From top to bottom: Average Lorentz Fig. 9a. Run 7: Spectrum for initial Lorentz factor of the 

factor per particle, for number density of the bunch n = beam particles 7 = 1000 (solid), 500 (dottedO), 250 (short- 

0.1, 0.2, 0.3, 0.4 «gj* (labelling as in Fig. I8a|) . Clearly, dashed), 100 (long-dashed). The central frequency shift as 

a small beam number density results in negligible losses, j 2 . 
and therefore, no coherent emission. 
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Fig. 9b. Pulse power for Run 7 plotted. 
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Fig. 9c. Run 7: The average Lorentz factor versus the av- 
erage distance z of the particles in simulations with chang- 
ing Lorentz factor. 



